Spin evolution of spin-1 Bose-Einstein condensates 
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An analytical formula is obtained to describe the evolution of the average populations of spin 
components of spin-1 atomic gases. The formula is derived from the exact time-dependent solution 
of the Hamiltonian Hs = cS^ without using approximation. Therefore it goes beyond the mean 
field theory and provides a general, accurate, and complete description for the whole process of non- 
dissipative evolution starting from various initial states. The numerical results directly given by the 
formula coincide qualitatively well with existing experimental data, and also with other theoretical 
results from solving dynamic differential equations. For some special cases of initial state, instead 
of undergoing strong oscillation as found previously, the evolution is found to go on very steadily in 
a very long duration. 
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The liberation of the freedoms of spin of atoms in opti- 
cal traps opens a new field, namely spin dy- 
namics of condensates, which is promising for super-high 
precise measurement, quantum cornputation, and quan- 
tum information processing. [1, 0, [1] Recently, the evo- 
lution of spinor condensates has been extensi vely studied 
experimentally and theoretically. [l^, [H, 113, Hi] Ini- 
tially, the condensate was prepared in a Fock-state or a 
coherent state confined in an optical trap. Then, due 
to the spin-dependent interaction, the system begin to 
evolve where a pair of atoms with spin components 1 
and -1 can jump to an 0, and vice versa, via scatter- 
ing. Finally the system will arrive in equilibrium, how- 
ever the process is not smooth. In 1998, the average 
population of each of the spin components fi =1, 0, and 
-1 was found to depend sensitively on initial states and 
may oscillate strongly with time. [l^ . This finding was 
further confirmed by a number of research groups. In 
2006, in the study of the probability of finding a given 
number of bosons in a given ^ state, the " quantum car- 
pet" spin-time structure was found. |l3| These findings 
show the amazing peculiarity of the spin dynamics. Re- 
lated theoretic calculations are mostly based on the mean 
field theory. Although, in a number of particular cases, 
theoretical results compares qualitatively well with ex- 
perimental data, the underlying physics remains to be 
further clarified. This paper is a study of the evolution 
of the average populations. We shall go beyond the mean 
field theory but use strict quantum mechanic many-body 
theory with a full consideration of symmetry. Instead 
of solving dynamic differential equations under specified 
initial condition, we succeed to derive a general analyt- 
ical formula to describe rigorously the whole process of 
evolution (non-dissipative) and is valid for all possible 
initial status. This is reported as follows. 

It is first assumed that the initial state of N spin-1 



atoms is a Fock-state with populations A^i, A'o and A^_i, 
the magnetization M — Ni — N^i. When N and M are 
given, the Fock-state can be simply denoted as |A''o). Let 
the part of the Hamiltonian responsible for spin evolution 
be Hs ~ cS^, where c is a constant, S is the operator of 
total spin. Then, the time evolution reads 
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where t = ct/h, and j^.j > is the all-symmetric to- 
tal spin-state with good quantum numbers S and AI. 
By using the analytical forms of the fractional parent- 
age coefficients and Clebesh-Gordan coefficients fll, [13, 
particle 1 can be extracted from the total spin-state as 



(2) 



where is the spin-state of particle 1. The coeffi- 

cients involved in ([1]) and Q are given in the appendix. 
Inserting ^ into ([T]), the probability of particle 1 in /i 
can be obtained, it reads 



(3) 



where 



^N.,, {No\^s,m){^Im\No) (4) 



P„^'^^ = {A{N, S, M, + {B{N, S, M, fi))' (5) 
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The summation covers S = N, N — 2, M* , where 

M* = M {oT M + 1) \i N - M is even (or odd). 
Since the particles are identical, each of them plays the 
same role, therefore the average population in fi is just 



M 
N„. 



j(r) = (a+a^) (this identity has been exactly 



NF 

proved numerically). In what follows fi — is assumed 
(the cases with /i / can be thereby understood). The 
label /i may be neglected from now on if /i = 0. 

Eq.Q is an exact consequence of the Hamiltonian 
Hg = cS^, no approximation has been introduced, it 
gives an analytical description of the whole evolution 
(non-dissipative) . There are time dependent and in- 
dependent terms, it implies an oscillation surrounding 
a background. It is straight forward from ^ that 
P^^{t) = P|fJ-T) = Pifjr + ^), therefore PjV^Jf + 
r) = Pji^ (^ — T).It implies that the oscillation is periodic 
with the period tt and Pj^ (r) is symmetric with respect 
to T = f . Furthermore, since cos(4(5' + 3/2)(f + r)) = 



- cos(4(5 + 3/2) (f - r)), 0^^{t) is antisymmetric with 
respect to f , we have Pj^Jf+r) = 2B]^^^ -P]lf^(f -r). 
Therefore, once Pj!f (r) has been known in the domain 
to 7r/4, it can be known everywhere. In particu- 
lar, PjV^jO) = 7Vo/iV, 
2B*f^ - iVo/iV. 
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In (4) the factor Pq'^^ has an exact analytical form as 



- M (2 + l/N)S{S + 1) - 1 - AP{2 + 3/N) 



{2S + 3)(2S-1) 



(8) 



When N is large, Pq 



^ 1 



(1 - (M/5)2). Therefore, 



In particular, when M — > 0, Bj!/^ « i. The value 1/2 
was first obtained numerically by Law, et al (l^ . and was 
supported by the recent study Ijy Chang, et al @. Now 
this value is obtained analytically, and is further found 
not depending on Nq. When M ^ N, S must also tend 



and BfJ^ 



as it should be. 



to N, therefore both Pf 

For the time-dependent term, O'^'^ in ^ depends on 
Nq strongly. There are three representative cases. 

(i) When Nq ^ N - M or 0, O^^f is distributed in 
a narrow domain of 5* (say, from Sa to Sb) as shown 
in Figd^ and [Tt. In this case, when O^'^ is roughly 
considered as a constant in the narrow domain, from ^ 
we have 



Oj^l(r) « J2 cos(4(2fc + Sa + 3/2)r) = P^'Gir) 

(10) 

where Pf^ is time-independent, k — (S ~ Sa)/2, fc,„ax — 
{Sb — Sa)/2. G{t) can be exactly rewritten as 



G(r) = cos(4(S'a+3/2+fc„,ax)r) sin(4(A:„iax+l)T)/ sin(4T) 

(11) 



The denominator sin(4r) affects the behavior of G(t) 
strongly. In the neighborhoods of 0, the magnitude of 
G(r) would be remarkably larger because sin(4T) is small, 
in particular, G(0) = fcmax -|- 1. In the neighborhoods of 
7r/4, the magnitude of G(t) would also be larger due 
to the denominator. However, since G(7r/4) = 0, there 
would be a strong oscillation when r — > 7r/4. 

(ii) When Na~ (N - M)/2, Of^f is distributed in a 
broad domain of S as shown in Fig[T}; where O^'^ and 



O 



M,S+2 
No 



have similar magnitudes but opposite signs. In 



this case, the summation in ([5]) can be divided into two, 
similarly we can define 



G{t) = ^ cos(4(4fc' + 5a + 3/2)T) 



k'=0 



- J2 cos(4(4fc" + Sa + 7/2)r) 

fc"=0 

2 sin(4T) 
sin(8T) ' 



(12) 



sin(4(S'a + 2 + 

max 

)t) Sin(8(fcniax + 1)^) 



The feature of G(t) is greatly different from G(t), in 

particular G(0) = G(7r/4) = 0, the denominator sin(8T) 

implies that G(t) would be large in the neighborhood 
of T « tt/S. This leads to a very different feature of 
evolution as shown later. 

(iii) When Nq is not close to the above cases, the vari- 
ation of Ojy ' against S has a band structure as shown in 

FigUJi, where neighboring O^'^ and O^'^^"^ may have 
the same or opposite signs. 

Examples of P]lf (r) calculated from ([3]) are given in 
the follows. Fig[2] shows the evolution in the whole pe- 
riod to TT, where the strong oscillation is concentrated 
in the neighborhoods of kn/A (a) or fc7r/4 -I- tt/S (b), 
k is an integer, due to the distinct features of G(r) 

and G{t). These figures show the symmetry in the 
period. Experimentally, the duration of observation is 
much shorter than tt. Evaluate under the Thomas-Fermi 
limit, when the trap is described by an isotropic har- 
monic potential with frequency a;/27r, r = tt is associated 
with tperiod = TT{N/uj^ f/^Xsec, where X = 1.52 x 10"* 
(3.86 xlO^) for ^^Rb (^^Na). In what follows r is only 
given in a short duration. 

The cases A^o — N — M are shown in Fig[3^ to[2l3. Fig. 
[3]i is associated with the experiments by the MIT group 
(upper panel of Fig. 2 of 3); Fig. [3Jd and c are the cases 
that experiment error emerges which makes M deviate 
from slightly. Fig. and e are associated with the 
experiments by GIT group (Fig.l of [l3| ), and Hamburg 
group (Fig. 5 of 0), respectively. Where, all Pj^ (r) (in 



solid lines) tend to B 
as predicted above. 



1/2 or lower (if M is larger) 
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The cases Nq = are shown in FiglSf to [3h, respec- 
tively. Where [3f is associated with the lower panel of 
Fig.2 of ref. 3 [Stenger98]. 

The cases Nq = {N - M) /2 are shown in FigO and 
[Sj. When M is small the evolution is very steady in a 
very long period to ^ 7r/8, then a strong oscillation 
occurs suddenly in the neighborhood of 7r/8 arising from 

the feature of G(t). Afterwards, the evolution becomes 
steady again, and repeatedly. 

When Nq is not close to the above cases, two examples 
are given in Figl51<; and Fig[3l. The former one is the 
case discussed by Law, et al. (shown in Fig. 3 of [[l^]). 
In this case, O^'^ is nearly chaos (FigUJi), (r) oscil- 
lates with T with a very high frequency in the beginning, 
but suddenly disappears, and suddenly recovers, and re- 
peatedly. 

In summary, this paper has essentially two findings 

(1) Going beyond the mean field theory, without the 
necessity to solve dynamical equations, a general ana- 
lytical formula has been derived based on symmetry to 
describe the evolution of the average populations P]!f (r) 
initiated from a pure Fock-state. This formula is an ex- 
act consequence of the Hamiltonian Hs — cS^ with a full 
consideration of symmetry, no approximation is adopted. 
Therefore the analysis based on this formula can help us 
to understand better the peculiarity of spin evolution. 
For examples, one can understand why the oscillation 
of P]!/ (r) becomes very strong in somewhere (in 7r/4 or 
7r/8), why (r) is symmetric with respect to 7r/2, and 
so on. The results from the formula coincides qualita- 
tively with existing experimental data or other theoreti- 
cal results. It is expected that, when accurate experimen- 
tal data come out, a detailed quantitative comparison can 
be made. 

(2) A special initial state with Nq = {N - M)/2 and 
Af « was found where the evolution of P]^ (r) is steady 
in a very long duration from the begining until r « 7r/8 
This special stability is noticeable. 

When the initial state is not a pure Fock-state but 
a superposition of them, the generalization is straight 



forward. 
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Appendix 

The coefficients in dH) and ^ are given as follows 

A{N, S, M, ^i) = af^ Cf^*S+i3/-M (13) 

B{N, S, M, ^i) = 6f C^^%^,^M-, (14) 

where 

4^' = [(1 + (-l)^"^)(iV - S)iS + 1)/{2N{2S + 

(15) 

C = [(1 + (-1)^^^) S{N + S+ l)/(2iV(25 + l))]i/2 

(16) 

and C'f^*5_(_jjy,^_^ are the Clebesh-Gorden coefficients, 
their analytical forms are given in [l^ . 

The set of coefficients {^smI-^o) are obtained by diag- 
onalizing the matrix of operator 

{Nl\S^\Nb) ^ AqSn^^.n, + A+Sn;^,n,-2 + A-Sn;^,n,+2 

(17) 

where Aq = + N + Nb + 2NNb - 2N^, A+ = 
y^Nbim - l)iN + M -Nb + 2){N -Nb + 2) and 
A_ = ^{Nb + l)(^h + 2){N + M- Nb){N -M- Nb). 



References 



[1] Tin-Lun Ho, Phys. Rev. Lett., 81, 742 (1998) 

[2] D.M. Stamper-Kurn, M.R. Andrews, A. P. Chikkatur, S. 
Inouye, H.-J. Miesner, J. Stenger, and W. Ketterle, Phys. 
Rev. Lett., 80, 2027(1998) 

[3] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H. -J. Mies- 
ner, A. P. Chikkatur, and W. Ketterle, Nature (London) 
396, 345 (1998). 

[4] A. Goelitz, T. L. Gustavson, A. E. Leanhardt, R. Low, 
A. P. Chikkatur, S. Gupta, S. Inouye, D. E. Pritchard, 
and W. Ketterle, Phys. Rev. Lett., 90, 090401 (2003) 

[5] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett., 94, 160401(2005) 

[6] A. Sorensen, L.-M. Duan, J.I. Cirac, and P. Zoller, Na- 
ture 409, 63 (2001) 

[7] L.-M. Duan, J.I. Cirac, and P. Zoller, Phys. Rev. A 65, 



033619(2002) 

[8] K. Molmer and P. Zoller, Phys. Rev. A50, 67(1994) 
[9] M.-S. Chang, CD. Hamley, M.D. Barrett, J.A. Saner, 
K.M. Fortier, W.Zhang, L. You, and M.S. Chapman, 
Phys. Rev. Lett., 92, 140403(2004) 
[10] M.S. Chang, Q. Qin, W.X. Zhang, L. You, and M.S. 

Chapman, Nature Physics (London) 1, 111(2005) 
[11] H. Schmaljohann, M. Erhard, J. Kronjager, K. Seng- 
stock, K. Bongs, Appl. Phys. B-Lasers and Optics 79, 
1001 (2004) 

[12] C.K. Law, H. Pu, and N.P. Bigelow, Phys. Rev. Lett., 
81, 5257(1998) 

[13] H. Pu, C.K. Law, S. Raghavan, J.H. Eberly, and N.P. 

Bigelow, Rhys. Rev. A., 60, 1463(1999) 
[14] R.B. Diener and T.L. Ho, arxiv:cond-mat /0608732 



4 



(2006) 

[15] C.G. Bao and Z.B. Li, Phys. Rev. A 70, 043620 (2004) 
[16] C.G. Bao and Z.B. Li, Phys. Rev. A 72, 043614 (2005) 
[17] C.G. Bao, Frontier in Physics in China 1, 92 (2006) 
[18] Z.B. Li and C.G. Bao, preprint, cond-mat/0703737l 
[19] A.R. Edmonds, Angular Momentum in Quantum Me- 
chanics (Princeton, New Jersey, 1960) 



5 



-0.0006 



10 20 30 

A g\ (a)N„=N-M 





400 



(d) N =(N-M)/3, M=0 



200 



400 600 

S 



800 



1000 



FIG. 1: O^f versus S. N = 1000 is given (the same in the 
follows) . 
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(a) N =N-M, M=0 
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FIG. 2: Evolution of the average population with n = 0. 
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FIG. 3: Evolution of the average populations with fi — 
(black), 1 (red), and —1 (blue). For the case M = 0, the red 
and blue lines overlap. 



